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AUTOMATIC EATA PROCESSING FOR PHOTOGRAPHIC 


PHOTOMETRY IN SPECTROGRAPHIC ANALYSIS 
by William A. Gordon and Anne K. Gallagher 
Lewis Research Center 


SUMMARY 

Calculations required to convert photographic densities of spectral lines 
to the concentration-related intensity parameter in spectrographic analysis 
have been programmed for digital computer processing. 

The program produced an emulsion calibration curve for each data set and 
converted all percent transmission values to relative intensities. Further, 
the line intensity ratios of the analytical lines to the internal standard 
lines were printed out for each step and exposure. These calculations were 
made for a maximum of 50 elements, four internal standard lines, and 11 expo- 
sures made with a seven-step sector. 

The program is written in Fortran IV (version 13) and is operational on 
the IBM 7094-2/7044 direct -coupled system of the Lewis Research Center. 

INTRODUCTION 

In emission spectrographic analysis the concentrations of the various ele- 
ments in a sample are determined from the photographic density of selected 
spectral lines that are recorded on a photographic plate. The conversion of 
the data from photographic density to concentration requires lengthy calcula- 
tions and graphical procedures that are adaptable to automatic processing. The 
description of the graphical construction of the emulsion calibration curve from 
step-sectored spectra follows a standard practice as described in reference 1. 
The processing of the most laborious steps in the calculations by an electronic 
computer is the subject of this report. 

The program, described in appendix B determines the emulsion calibration 
curve for each photographic plate and converts all percent transmissions to the 
intensity parameters from which element concentrations are derived. Calcula- 
tions were made for a maximum of 50 elements, in a maximum of 11 samples. In 
spectral analysis using photographic photometry it is customary to stepwise 
attentuate the intensity of the spectral lines along their length to calibrate 
the emulsion and also to extend the concentration range for each exposure. 

This program allows the recording of data in a maximum of seven steps of atten- 






uation for each spectral line. Conventionally, the element line intensities 
are ratioed to an internal standard line intensity to compensate for some ex- 
perimental errors in the procedure. In this program provisions were made for 
calculating these ratios using a maximum of four internal standard line inten- 
sities. This was done to facilitate selection of internal standard lines in 
the development of analytical methods. A typical set of calculations, requir- 
ing about 40 hours when performed manually, may be completed in about 1 minute 
of computer time. In addition, the automatic calculations are less subject to 
error than those performed manually. 


PHOTOGRAPHIC DATA REDUCTION 

The conversion of photographic density to element concentration can be ar- 
bitrarily divided into four basic operations. Two of the operations are per- 
formed manually and two by computer as indicated: 

(1) The conversion of photographic density into percent transmission using 

a densitometer (Manual) 

(2) Graphical determination of the relation between percent transmission 

and the spectral line intensity that produced the photographic den- 
sity, referred to as the emulsion calibration curve (ECC) (Computer) 

(3) Calculation of the ratios of the analytical line intensities and the 

internal standard intensity, with appropriate background corrections 
(Computer) 

(4) Graphical conversion of the line intensity ratios to elemental concen- 

trations by means of calibration curves prepared from standards for 
each element sought (Manual) 

The automation of the second operation, which involves the plotting of 
experimental data followed by two separate curve -fitting procedures, provides 
the emulsion calibration curve ECC. Although the mathematical relation for an 
emulsion calibration curve can be approximated in general by a mathematical 
function, a more precise procedure is to determine experimentally the ECC for 
each photographic plate. It should also be noted that the ECC is a function 
of the wavelength of the exciting light and for highest precision should there- 
fore be determined for selected wavelength intervals. In the present work, 
however, a single ECC is determined representing an average over the experi- 
mental wavelength region. Additions can be made to this program to take into 
account this wavelength dependency should this be desired. 

In the first curve -fitting procedure in operation ( 2 ), the percent trans- 
mission readings of adjacent steps of each spectral line recorded on card 3 
(appendix A, fig. 1(b)), are input on linear X-Y coordinates. These ordered 
pairs (X,Y) have a definite intensity relation, known as the "step factor," 
which is determined by the physical configuration of the step sector used in 
the procedure to stepwise attenuate the light. In this work the step factor 
is 2.0 (i.e., any two adjacent steps of a given line have an intensity rela- 
tion of 2.0). Similarly, alternate steps for a given line have an intensity 
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ratio of 2 , and for steps separated by n step positions the intensity ratio 
of the darker step to the lighter step is 2 n . Thus, for the preliminary curve 
each ordered pair (X,Y) is derived from transmission readings for n * 1 (i.e., 
adjacent steps in the same line). It should be noted that these ordered pairs 
(X,Y) are derived from line, not background readings. 

The program, as designed, orders a maximum of 250 pairs (X,Y) from the in- 
put data including two pairs fixed at (0,0) and (100,100), which are weighted 
to force the fit. Should the maximum number of ordered pairs exceed 250, it 
may be necessary to revise card 1 (appendix A, fig. l(a)), so that the ordering 
of the elements will provide representative data with respect to wavelength re- 
gion and to the range of densities in the first 250 ordered pairs. After cali- 
bration all of the input data are used in the ratio calculations. 

The total number of points registered are curve fitted to a polynomial of 
the form 


Y = A 1 + A 2 X + A 3 X 2 + A 4 X 3 

by a general least-squares method using the Crout reduction and a method devel- 
oped at the Lewis Research Center by Bert Henry to invert the matrix. The re- 
sulting coefficients describe a curve known as a "preliminary emulsion calibra- 
tion curve." This curve-fitting operation serves to effectively smooth all the 
plotted data. 

The second curve -fitting operation results in the final emulsion calibra- 
tion curve ECC, which relates percent transmission (ordinate) to the intensity 
parameter (abscissa). Values for plotting this curve are derived from the pre- 
liminary curve as follows: start at X = 0.02 (2 percent transmission) on the 

preliminary curve, and compute Yp; now let Yp = Xg, and find Yg correspond- 
ing to this Xg from the preliminary curve function. Then set Yg = X3, and 
compute Y3, • • • , etc. In this way a series of X's (maximum of 10) is gen- 
erated, which represent the smoothed percent-transmission readings whose inten- 
sity relation are known from the step factor as explained previously. In the 
experiments presented herein this relation may be expressed as the percent 
transmission T = 2 n Z, where n is the difference in step number and Z is 
the intensity parameter that is related to element concentration in analytical 
work. The conversion of all percent -transmission data to this intensity param- 
eter is done by a Lagrange four-point interpolation technique. The interpola- 
tion was easily performed on a large-scale computer and yielded greater accura- 
cy than the method of least squares. 

The programming of operation (3) requires the straightforward calculations 
of spectral line intensity ratios using the intensity parameters derived from 
the percent transmissions and the ECC, as described in operation (2). All in- 
tensity ratios were tabulated in both the background-corrected and the non- 
corrected forms; however, when background readings were not entered on card 3 
(appendix A, fig. 1(b)), for either the internal standard or the analytical 
line, a zero was printed in the column reserved for the background-corrected 
ratio, as shown in the sample output for nickel (appendix C). In the calcula- 
tion of line intensity ratios instructions are programmed so that maximum use 
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is made of the available data. For example, percent -transmission readings in a 
given exposure may not be recorded in the same step for the internal standard 
and the analytical lines. When this occurs, a calculation is made that auto- 
matically corrects for any step difference between the internal standard and 
the analytical lines. This calculation is performed by using 2 n as a multi- 
plier, where 2 is the known step factor and n is the step difference. The 
intensity ratios obtained by this procedure may be subject to extrapolative er- 
rors; therefore, it is desirable to identify such values to warn of possible 
error by printing an asterisk to the right of the intensity ratios obtained by 
extrapolation. It should also be noted on the sample output (appendix C) that 
an average value is printed for each element and exposure. In the calculation 
of this average the extrapolated intensity ratios are given the same weight as 
the other intensity ratios. Although this may not be a statistically valid 
procedure, the error due to this simplification is not significant for most 
work at the Lewis Research Center. 

Since a considerable amount of manual effort is still necessary to record 
the input data, errors in input can be expected to arise with about the same 
frequency as when calculations are made manually. The error-searching proce- 
dure is performed manually and is similar to the general procedure used by the 
spectrographer when making manual calculations. The preliminary curve is 
printed out, and wide deviations from the average are often recognized as erro- 
neous data. The dubious value can then be traced back through the input, 
reread on the densitometer, and corrected if necessary. In addition, the per- 
centage deviation of each of the data points from the preliminary curve fit is 
printed out in tabular form and serves as a useful means for revealing ques- 
tionable data. 

The line intensity ratios calculated by the program are then converted 
manually to element concentrations by means of calibration curves prepared from 
standards for each element sought. These calibrations cannot be conveniently 
programmed because the parameters under study at the Lewis Research Center af- 
fect the positions of the concentration curves. Once these calibrations are 
permanently established, the last operation, conversion of line intensity ra- 
tios to concentrations, can be included in the computer program by using curve 
fitting techniques similar to those described previously. 

This program for spectrographic calculations forms a basis to which addi- 
tional modifications can be made to make it more comprehensive. In its present 
state , the program serves the needs of a spectrographic laboratory engaged in 
both research and service analyses. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 23, 1965. 
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APPENDIX A 


INPUT 

This data reduction program is written in Fortran IV and is operational on 
the IBM 7094-2/7044 direct coupled system of the Lewis Research Center. The 
main program and its eight subroutines require 8000)g store in common data, 
30000 ) 8 store for the program, and 25000)g additional store for Lewis monitor 
systems routines. For a smaller computer installation, overlay or some other 
store-sharing technique would he necessary. 

Figure 1 contains a list of the input with samples of the forms used at 
Lewis. The input in card type 3 is sorted in columns 11 to 7 to provide for 
the proper pairing of data in the preliminary • curve fit. 

The sample data (fig. l) and the sample output (appendix C) are abbrevi- 
ated to show only the first two internal standards and the element nickel. 

Instructions for preparing the input are given in table I. 


TABLE I. - INSTRUCTIONS FOR PREPARING INPUT 
[See listing for specific formats.] 


Card, type 

Number 

Columns 

Description of variables 

1 

( Element 
data) 

1 

NEL 

1 to 5 

7 to 12 
19 to 24 

NEL number of elements present 

(maximum, 50 ) 

EL 6 letter code for element 

ELB 6 letter code for element 

background 

2 

(Constants for 
one case) 

1 

2 to 5 

10 to 15 
16 and 17 

18 and 19 

NMAX number percent transmissions 

recorded (maximum, 3850) 
HATE plate number 

KSTOP number standards to be used 

(maximum, 4) 

NEXP number exposures (maximum, ll) 

1 

1 to 72 

(DATE (i), any statement about the data 
1=1,12) 

3 

(Data for 
one case) 

NMAX 

7 and 8 

9 and 10 
11 

13 to 18 
19 to 24 

I element number (must match order 

of type 1, 1 to 4 standards) 

J exposure number 

K step number (maximum, l) 

CLU percent transmission 

BAK percent transmission of 

background 


Card types 1 and 3 may be repeated for a number of cases, but card 
type 2 must be the same for all cases in a given run. Forms used at Lewis 
for card types 1 and 3 are shown in figure 1. 
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(a) Card type 1. Background measurements: L is measured (b) Card type 3. 

on low wavelength side of line; H is measured on high 
wavelength side of line; A is measured as average of both 
sides of line. 

Figure 1. - Sample input 







APPENDIX B 


FORTRAN IV PROGRAM 

A synopsis of the subroutines and their content listed in order of calling 
sequence is as follows: 

BEGIN MAIN PROGRAM 

(1) Inputs data 

(2) Sets up X,Y-pairs for PC curve fit 

(3) Outputs relative intensities 

GLSQl Performs a least-squares polynomial fit on X,Y-pairs 

PLOTXY Plots preliminary curve ; PLOTXY is Lewis Monitor Subroutine ( see 
ref. 2) 

FUNC Function Subprogram to generate independent variable function 

EFUNC Function Subprogram to backtransf orm independent variable function 

DINVER Subroutine using Crout reduction and a routine to invert matrix 

DMVMUL Subroutine to obtain solution vector from the inverse 

CALC Subroutine to determine X,Y-pairs for the ECC from the fitted PC 

LGRNGE Subroutine to perform four-point Lagrange interpolation; converts per- 
cent transmission to relative intensity 

RATIO Subroutine which ratios relative intensity of an element to that of a 
standard and corrects for background if necessary (maximum number of 
standards, 4) 

Flow charts of BEGIN, CALC, and RATIO are shown in figure 2. These three 
routines are pertinent to the spectrographic data analysis. The remaining sub- 
routines are in general use at the Lewis Research Center. The main program and 
subroutines used for data reduction are presented after figure 2. 
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$ 1 BF TC BEGIN LIST.OECK 

AUTOMATIC DATA PROCESSING FOR SPECTROGRAPH IC CALCULATIONS 
COMMON I W, KSTOP , NEXP 
COMMON Z(50,ll,7),ZB(50,ll,7) 

COMMON EL!50),£LB(50) 

DIMENSION VAL! 50, II, 71 ,8AKVAL<50, 11, 7) ,X(250) ,Y(500),P(35) 
DIMENSION ZX(200) ,ZY< 200) 

DIMENSION DATE! 12) 

EQUIVALENCE (Z.VAL) , ( ZB» BAKVAL ) 

READ (5*98) NEL 
DO l 1=1, NEL 

1 READ! 5, 102 ) EL( I),ELB(I) 

1000 WRITE! 6, 106) 

READ! 5, 100) NMAX, PLATE, KSTOP, NEXP 
RE AD! 5 , 99 ) t DATE! I ), 1=1, 12) 


NMAX= NUMBER READINGS ON THIS PLATE 

NEL= NUMBER OF ELEMENTS PRESENT IN THIS SAMPLE 
ELEMENTS INPUT FIRST (MAX=4) ARE TO BE USED AS INTERNAL STANDARDS. 
AG AND AGX = SILVER, A=ARGUN, W = TUNGSTEN. 


EL AND ELB ALLOW A MAXIMUM 6 CHARACTER DESIGNATION FUR EACH 
ELEMENT AND ITS BACKGROUND TO BE USED AS COLUMN HEADINGS. 
VALUE FOR - I - MUST MATCH ITS CODES IN EL! I), ETC) 

IF (NMAX— 250)3,3,2 

2 NSTOP = 250 
GO TO 4 

3 NSTOP = NMAX 

4 JJ=0 
III = 0 
KK = 0 
NOM = 1 

BEGIN DETERMINATION OF 
DO 5 LK=1, NSTOP 


POINTS TO BE USED IN FINDING THE PC 


BEGOOOOO 
BEG00010 
BEG00020 
8EG00D30 
BEG00340 
BEG00050 
BEG00060 
BEGOOO 70 
BEG00D80 
BEG00D90 
BEG00100 
BEG00110 
BEG001 20 
BEG00130 
BEG00140 
BEG00150 
BEG00160 
BEG00170 
BEG00180 
BEG00190 
BEG00200 
BEG00210 
BEG 220 
BEG 230 
BEG 240 
BEG00250 
BEG00260 
BEG002 70 
READIN BEG00280 
BEG00290 
BEG00300 
BEG00310 
BEG00320 
BEG00330 
BEG00340 
8EG00350 
BEG00360 
BEG00370 
BEG00380 
8EG00390 
BEG00400 


DO 10 I I = NOM , NMAX 







BEG004 10 
BEG00420 

1= ELEMENT NUMBER, XX , 

INTEGER 

VALUE 

FROM 

l 

TO 

50 

BEG00430 

J=EXPUSURE NUMBER, XX, 

INTEGER 

VALUE 

FROM 

1 

TO 

11 

BEG00440 

K= STEP NUMBER, X, 

INTEGER 

VALUE 

FROM 

1 

TO 

7 

BEG00450 


CLU = PERCENT T 

BAK = PERCENT T DUE TO BACKGROUND. 

READ (5,101) I , J , K , CLU, BAK 
VAL! I, J,K) = CLU 
BAKVAL! I,J,K) = BAK 
IF (JJ-J) 8,6,8 

6 IF (KK-K) 8,7,8 

7 Y( LK) = VAL! I, J,K) 

X ( LK ) = VAL! I, J,K-l) 

DEBUG LK , X ( LK ) , Y ( LK ) 

NOM = II+l 

JJ = J 
KK = K+i 
II I = I 


BEG00460 
BEG00470 
BEG00480 
BEG00490 
8EG00500 
BEG005 10 
BEG00520 
BEG00530 
BEG00540 
BEG00550 
BEG00560 
BEG00570 
BEG00580 
8EG00590 
BEG00600 
BEG00610 
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non 


GO TO 9 

8 1 1 1= I 
JJ = J 
KK = K + l 
N0M = 1 1 +1 

IF ( NOM-NMAX ) 10* 10,9 

10 CONTINUE 

9 N0=LK 

IF (NOM-NMAX) 5,11,11 
5 CONTINUE 

11 CONTINUE 

C POINTS ARE AVAILABLE. CALL GLSQ1 TO DETERMINE PC 
IW=l 

CALL GL SQ1 ( X, Y, NO, P ) 

C PC HAS BEEN DETERMINED. ECC HAS BEEN DEFINED. 

WR ITE ( 6, 106) 

WRITE! 6, 104) PLATE 
WRITE! 6, 105) 

WR I TE ( 6 , 107 ) (DATE! I ), 1 = 1, 12) 

KODE= 1 
C 

C KODE = l INDICATES DERIVATIVE IS NOT TO BE CALCUALTED IN SUB. LGRNGE 
C RENAME RESULTS OF SUB. CALC 
DO 15 J = 1 , NU 
ZX! J)=X! J) 

15 ZY( J ) = Y( J ) 

C TRANSFORM PERCENT TS TO RELATIVE INTENSITIES BY MEANS OF LAGRANGE 
C INTERPOLATION. 

DO 30 J= 1 » NEXP 
DO 25 K= 1 , 7 
DO 20 1=1, NEL 

IF ( VAL! I, J,K)-1. ) 16, 16, 17 

16 Z( I, J,K) =0. 

ZB! I, J,K)=0. 

GO TO 20 

17 CALL LGRNGE(ZX,ZY,NO,VAL( I , J , K ) , Z ( I , J , X ) ,DER,KODE) 

IF ( BAKVAL! I, J,K)-1. ) 18, 18, 19 

18 ZB! I , J , K ) =0. 

GO TO 20 

19 CALL LGRNGE (ZX,ZY,NO,BAKVAL(I,J,K),ZB(I,J,K),DER,KODE) 

20 CONTINUE 
25 CONTINUE 

30 CONTINUE 
L= 1 
LL=6 

PRINT RELATIVE INTENSITIES 

31 WR ITE ( 6, 108 ) ( EL ! I ) , ELB! I ) , I =L , LL ) 

DO 50 J=1 , NEXP 
WRITE! 6, 103) J 
DO AO K=1 , 7 

WRITE! 6,109) K, ( Z ( I , J , K ) , ZB ( I , J , K ) , I =L , LL ) 

AO CONTINUE 

50 CONTINUE 
LL=LL-*-6 
L=L»6 

IF ( LL— NEL )51,51,52 

51 GO TO 31 

C BEGIN FINDING RATIOS OF RELATIVE INTENSITIES. 

52 CALL RATIO (NEL, Z, ZB) 

WR ITE! 6, 106) 


8EG00620 

BEG00630 

BEG006A0 

BEG00650 

BE G00660 

BEG00670 

BEG306B0 

8EG00690 

BEG00700 

3EG00710 

8EG00720 

BEG00730 

BEG007A0 

BEG00750 

BEG00760 

BEG00770 

BEG00780 

BEG00790 

BEG00800 

BEG00810 

BEG00820 

BEG00830 

8EG008A0 

BEG00850 

BEG00863 

BEG00870 

BEG00880 

BEG00890 

BEG00900 

BEG00910 

BEG00920 

BEG00930 

BEG009A0 

BEG00950 

BEG00960 

BEG00970 

BEG00980 

BEG00990 

BEG01000 

BEG01310 

BEG01020 

BEG01030 

BEG010A0 

BEG01050 

BEG01060 

BEG01070 

BEG01080 

BEG01090 

BEG01100 

BEG01110 

BEG31120 

BEG01130 

BEG011A0 

BEG01150 

BEG01160 

BEG01170 

BEG01180 

BEG01190 

BEG01200 

BEG01210 

BEG01220 

BEG01230 
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non 


00 70 1=1,7 

00 65 J=l,ll 

00 60 K=1 , 50 

VALU.J, I)=0. 
BAKVAL ( K, J , I ) =0. 
60 CONTINUE 
65 CONTINUE 
70 CONTINUE 
GO TO 1000 


98 

99 
100 
101 
102 

103 

104 


FORMAT STATEMENTS 


, 2 ) 


FORMAT (151 
FORMAT C12A6) 

FORMAT (X,I4,4X,A6,2I2) 

FORMAT (6X, 212, U,X,2F6 
FORMAT (6X,A6,6X,A6) 

FORMAT ( 4H EXP, 12) 

FORMAT! 1H2, 35X,45H SPECTROGRAPH IC 
XA6 ) 

105 FORMAT! IH0,30X,62H PER CENT TRANSMISSION 
XEL INTENSITY ) 


CALCULATIONS FOR PLATE NO. 


READINGS CONVERTED TO 


8EG01240 
BEG01250 
BEGO 1260 
BEG01270 
BEG01280 
BEG01290 
BEG01300 
BEG01310 
8EG01320 
BEG01330 
BEGO 1340 
BEG01350 
BEG01360 
BEG01370 
BEG01380 
BEG01390 
BEG01400 
BEGO 1410 
, BEG01420 
BEG01430 
RBEG01440 
BEG01450 


106 

FORMAT! 1H1) 

BEG01460 

107 

FORMAT! 1H0.29X, 12A6) 

8EG01470 

108 

FORMAT! 1H1,6< A6,5X,A6,4X) ) 

BEG01480 

109 

FORMAT ( 1H , 12, ! 12G10.3) ) 

BEG01490 


STOP 

BEG01500 


END 

BEG015 10 
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$1 BFTC GLSWR L I ST , DECK, OEBUG 

SUBROUT INE GLSQ11X.Y, NSMALL, P) 

C 

COMMON IW 

DIMENSION XI 250) , Y( 500) , U( 351 ,A( 35 , 35 ) , ASA VE < 35 , 35 ) 

1. PI35J ,Pl(35) , YCALCI 250) .DELYIZ50) , XP ( 250) , YP I 500 > , ERRATA I 250 ) 

2, £LM(2) , TITLE! 1 2 ) , FMT < 1 2 ) , FMTOI 1 2 > , W! 250 ) 

DIMENSION AIN!35,35) 

DIMENSION FMTW(7) ,KMTW! 7) ,PK!1) ,KPL( 15) ,CC! 35) 

EQUIVALENCE ( FMTW! 5 I , KMTrt! 5 ) ) 

DATA (FMTWII), 1=1 , 7 ) /6H( 22H0 , 6HTHE P A , 6HRAMET E , 6HRS (A0.6H-A 12, 
16H5H) AR.3HE ) / , KMT 5 1 , KM T52/ 6H- A I1..6H-A 12,/ 

DATA(TITLE( I), 1 = 1,12) /6HAUT0MA, 6HT IC DA , 6HT A PRO , 6HCES S I N , 6HG FOR 
X,6HSPECTR, 6H0GRAPH* 6H I C C AL , 6HCUL A T I , 6HONS ,6H ,6H / 

MAIN PROGRAM FUR LEAST SQUARES ANALYSIS HAS BEEN ALTERED FOR THIS 
SPECIFIC PRUBLEM. 

ALL THE OPTIONS IN THE ORIGINAL PROGRAM HAVE BEEN REMOVED. 

DOUBLE PRECISION A , D , CC , UXMA X , A I N 
PK ( I ) = I .0 
KPL=64 
KPL! 2) =2 
POOL1=0. 

P00L2=0. 

NFMTQ=NBIG*J - 2 
NBIG=4 
I PL- 1 
I Y=0 

NBIGSV=N6IG 


NSMALL=NU. 

THE ORIGINAL 
IF IT BECOMES 
NNN=N8IG 
MAX =35 

DO 32 1=1, NSMALL 

32 W!l)=l. 
NSMALL=NSMALL+2 
NN=NSMALL 

X( NN-1 ) =0. 

Y( NN-1 ) =0. 

XI NN) = 100. 

Y( NN ) = 100. 

W! NN ) = 50 . 

W( NN-1 ) =50. 

KKK=1 

SN=NSMALL 

6000 DO 6010 1 = 1, NSMALL 
XP< I ) =X( I) 

6010 YP ( I)=Y( I) 

NBIG=NB IGSV 
NSAME=NSAMS V 
SUMYP=0 • 

XMAXX=0. 

DO 501 1=1, NSMALL 
IF ( XMAXX-A8S! X! I ) ) ) 

33 XMAXX=A8S( XI I ) ) 
DXMAX=XMAXX 

501 CONTINUE 


OF DATA POINTS NB I G=NO. OF COEFFICIENTS DESIRED 

DATA IS PLOTTED 

DESIRABLE TO OMIT THE PLOT CHANGE IPL TO D 


33,501,501 


GLSDOOOO 
GLSOOOIO 
GLS00D20 
GLS00030 
GLS00040 
GLS00050 
GLS00060 
GLS00070 
GLS00D80 
GLS00090 
, GLS00100 
GLS001 10 
GLS00120 
GLS00130 
GLS00140 
GLS001 50 
GLS00155 
GLS00160 
GLS00170 
GLS00180 
GLS00190 
GLS00200 
GLS002 10 
GLS00220 
GLS00230 
GLS00240 
GLS00250 
GLS00260 
GLS00270 
GLS00280 
GLSD0290 
GLS00300 
GLS003 10 
GLS00320 
GLS00330 
GLS00340 
GLS00350 
GLS00360 
GLS00370 
GLS00380 
GLS00390 
GLS00400 
GLS00410 
GLS00420 
GLS00430 
GLS00440 
GLS00450 
GLS00460 
GLS00470 
GLSQ0480 
GLS00490 
GLS00500 
GLS00510 
GLS00520 
GLS00530 
GLS00540 
GLS00550 
GLS00560 
GLS00570 
GLS00580 
GLS00590 
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OOO O O o ooo o o 


20 

1903 


22 

23 


35 

36 
2 


37 

38 


397 

300 


301 


130 


00 20 1 = 1* NSMALL 


GLS00600 

X< I »=X{ D/XMAXX 


GLS00610 

xp( i ) = x t n 


GLS00620 

yp( n=y< n 


GLS00630 

SUMYP=SUMYP+YP( I ) 


GLS00640 

SUMY=0. 


GLS00650 

CN = NB IG 


GLS00660 

KN=NBIG-1 


GLS00670 

DETNUM=0. 


GLS00680 

0ET0EN=0. 


GLS00690 

YPMEAN=SUMYP/SN 


GLS00700 

YMEAN= YPMEAN 


GLS007L0 

GLS00720 

GLS00730 

00 23 1=1, NSMALL 


GLS00740 

IF ( I Y 1 22,23,22 


GLS00750 

Y { I)=FUNC( I Y , Y ( I) ) 


GLS00760 

SUMY=SUMY+Y( I ) 


GLS00770 

YMEAN=SUMY/SN 


GLS00780 

GLS00790 

SET UP THE MATRIX OF 

COEFFICIENTS 

GLS00800 
GLS008 10 

DO 2 I=l,NBIG 


GLS00820 

DO 2 J = 1 » NB IG 


GLS00830 

A ( I, J)=0. 


GLSD0840 

00 2 K=l, NSMALL 


GLS00850 

F I =FUNC( l , X ( K ) I 


GLS00860 

FJ=FUNC ( J , X ( K ) ) 


GLSD0870 

IF (IWI 35,36,35 


GLS00880 

A ( I, J)=A< I, J)+FI*FJ*WCK) 


GLS00890 

GO TO 2 


GLS00900 

A( I, J)=A( I,J)+FI*FJ 


GLS00910 

CONTINUE 


GLS00920 

GLS00930 

DO 14 I=1,NBIG 


GLS00940 

D( I)=0. 


GLS00950 

00 14 K=l, NSMALL 


GLS00960 

FI=FUNC( I , X ( K ) ) 


GLS00970 

IF ( 1 W ) 37, 38,37 


GLS00980 

0( I)=D( I)*Y(K)*FI*W(K) 


GLS00990 

GO TO 14 


GLS01000 

0( I)=0( I ) +Y ( K ) *F I 


GLS01010 

CONT INUE 


GLS01020 

GLS01030 

WRITE (6,397) TITLE 


GLS01D40 

FORMATt 1HI, 10X, 12A6) 


GLS01050 

GLS01060 

CALL OINVER ( A , NNN , A I N , M AX 

, DET, ERR , D) 

GLS01370 

CALL DMVMUH A IN , NNN, 0, CC , MAX ) 

GLS01080 

XXAXX=1./XMAXX 


GLS01090 

DO 301 J=1,NBIG 


GLS01100 

P ( J ) =CC l J ) 


GLSOlllO 

P( J)=P( J)*FUNC( J.XXAXX) 


GLS01120 
GL SOI 1 30 

COMPUTE Y-CALC, OELTA-Y, AND 

STANDARO DEVIATION 

GL SO 1140 
GLS01150 

CONT INUE 


GLS01160 

DEBUG (P(J) , J=1,NBIG) 


GLSOl 170 

ERS=0. 


GL SO 1 1 80 

DEV=0. 


GL SO 1 1 90 

DO 8 1=1, NSMALL 


GLS01200 

X( I )=X( I )*XMAXX 


GLS012 ID 
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K=NSMALL+I GLS01220 

Y( K ) =0 . GLS01230 

DO 9 J=1,N8IG GLSO 12 AO 

FJ=FUNC ( J t X ( I ) ) GLS01250 

DEBUG FJ GLS01260 

9 Y(K)=Y(K)+P(J)*FJ GL SO 1 2 70 

TEM= Y( I ) -YMEAN GLS01280 

DETDEN=DETDEN«-TEM*TEM GLS01290 

TEM=YU)-YMtAN GLS01300 

DETNUM=DETNUM«-TEM*TEM GLS01310 

DELY( I > = Y C I > — Y ( K I GLS01320 

ERRATA( I)=DELY( I)/Y(K)*100. GLS01330 

ERS=ERS + ERRATA( I )»*2 GLS013AQ 

8 DEV=DEV+DELY( I )«DfcLY(I) GLS01350 

POOL 1=P00L 1 + DEV GLS01360 

ENDI V=SN-CN— l. GLS01370 

TVRDEN=TVRUEN+ENOIV GLS0I380 

TVRNUM=POUL 1 GLS01390 

100 DE V=DE V/ END I V GLS01A00 

DVTN=SQRT ( OEV) GLS01A10 

ERS=ERS/SN GLS01 A20 

ERA-SORT ( ERS) GLS01A30 

TQTVAR=TVRNUM/TVRDEN GLSQIAAO 

TOTDEV=SQRT( TOTVAR) GLS01A50 

DETRM=DETNUM/DETDEN GLS01A60 

CORRL = SQRT ( DETRM ) GLS01A70 

WR ITE( 6, AO 1 ) GLS01A80 

DO 10 1=1, NSMALL GLS01A90 

K=NSMALL«-I GLS01500 

10 WRITE! 6 t A02 ) X(I),YI1I,Y I K) , DELY! I ) , ERRATA ! I ) GLS01510 

MR ITE ( 6, A06 ) GLS01520 

A06 FORMAT ( A 1H0 THE REGRESSION EQUATION FOR THE ABOVE IS /IH0.35H Y = AGLS01530 

XO ♦ Al*X *■ A 2*X**2 * A3*X»*3 ) GLS015A0 

A09 IF (KN-10) A07,A05,A05 GLS01550 

AO 5 KMTWl 5 ) =KMT52 GLS01560 

GO TO A08 GLS01570 

A07 KMTWI 5 ) =KMT51 GLS01580 

A08 MR ITE (6,FMTW) KN GLS01590 

WRITE !6,A03> ( P ( l ) , I = 1 , NB IG ) GLS01600 

WRITE ( 6* AOA) DEV, DVTN , DETRM, CORRL , TOT VAR , TOTUEV, ERS, ERA GLS01610 

WRITE (6,91) OET GLS01620 

91 FORMAT ( 13H0DETERMINANT* GIA.6//) GLS01630 

IF(IW) A22 , A19 , A22 GLS016A0 

A 22 WR I TE( 6 , A30 ) GLS01650 

00 A20 1=1, NSMALL GLS01660 

IF (W(I)-l.) A2 1 , A20 , A21 GLS01670 

A21 WRITE(6,A31)X( I),Y(I),W( I) GLS01680 

A20 CONTINUE GLS01690 

A19 CONTINUE GLS01700 

A3 l F0RMATI1H ,2X,2HX=G15.6,2X,2HY=G15.6,2X, 12HWITH WEIGHT* F7.0) GLS0L710 

A30 FORMAT ( AAHOTHE ABOVE DATA HAS BEEN WEIGHTED AS FOLLOWS » GLS01720 

AO 1 FORMAT ( IHO,2X,10H IND. VAR. ,8X,10H DEP.VAR. ,8X,12HCALC. FUNC. GLS01730 

X,6X, 10H DEVIATION, 8X, 15H PER CENT ERROR ) GLS017AO 

AO 2 FORMAT ( 1H , ( 5( F12. A, 6X ) ) ) GLS01750 

AO 3 FORMAT ( 7GI 8 • 8 ) GLS01760 

AOA FORMAT! 1H0, 2X, 1AH THE VAR I ANCE=G1 5 . 7, 20H STANDARD DEVI ATI0N=G15. 7/GLS01770 
1 3X , 1AHDETERM INAT I0N=G15.7,8X, 12HC0RREL AT 1 0N=G1 5 . 7/ 1 7H POOLED VAR I AGLS01780 
2NCE=G15.7,3X, 17HP00LED OEV I A T I0N=G1 5. 7/ GLS01790 

3 1H0, 2X, 1AH PCT VARI ANCE=G15 . 7 , 20H STD. PCT DF VI AT 1 0N = G1 5 . 7 ) GLSO 1800 

GO TO ( 101, 102) ,KKK GLS01810 

101 DO 2A 1=1, NSMALL GLS01820 

K=NSMALL+I GLS01830 
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IF ( IY) 25,26,25 

GLS01840 

25 YP(KJ=BFUNC( IY, YIK) 1 

GLS01850 

GO TO 24 

GLS01860 

26 YP ( K ) =Y IK) 

GLS01870 

24 CONTINUE 

GLS01880 

IF ( IPL-I) 29,803,801 

GLS01890 

803 KKK=2 

GLS01900 

801 CONTINUE 

GLS01910 

KPLI 31 =NSMALL 

GLS01920 

K2N=NSMALL+NSMALL 

GLS01930 

102 WRITE(6,77)TITLE 

GLS01940 

77 FORMAT ( 2HPT, I2A6/2HPT/ 2HPT , 59X, 7HN0TE — /2HPT/2HPT, 59X,28H0R I GI NALGLSD1950 

1 DATA PLOTTED WITH • /2HPT, 59X, 30HCALCUL ATED 

DATA PLOTTED WITH •- /GLS01960 

X2HPT/2HPT, 59X, I4HPL0T OF PC ) 

GLS01970 

CALL SCALE ( NSMALL , X, KKLU ) 

GLS01980 

CALL SCALE ( K2N, Y , KKL Y ) 

GLS01990 

CALL PLOTMYIX ,Y .KPL.PK) 

GLS02000 

WR I TE( 6* 78 ) 

GLS02010 

78 FORMAT (2HPL,59X,23HLIGHT STEP VS DARK STEP 

) GLS02020 

802 CONTINUE 

GLS02030 

GO TO I 103,29) ,KKK 

GLS02040 

103 CONTINUE 

GLS02050 

27 DEV=0. 

GLS02060 

0ETNUM=0. 

GLS02070 

DETDEN=0. 

GLS02080 

DO 28 1=1, NSMALL 

GLS02090 

K=NSMALL +1 

GLS02100 

X( I)=XP( I) 

GLS02110 

Y ( I ) =YP ( I ) 

GLS02120 

Y(K)=YP(K) 

GLS02130 

TEM=Y( I ) -YPMEAN 

GLS02140 

DETDEN=DETDEN+TEM*TEM 

GLS02150 

TEM=Y(K)-YPMEAN 

GLS02160 

DETNUM=DETNUM+TEM*TEM 

GLS02170 

DELY ( 1 ) = Y ( I)-Y(K) 

GLS02180 

ERRATA ( I )=DELY( I ) / YP 1 K » 

GLS02190 

28 DEV=DEV+DELYU) *DELY(I) 

GLS02200 

POOL 2= POOL 2 ■••DEV 

GLS02210 

TVRNUM=P00L2 

GLS02220 

KKK=2 

GLS02230 

GO TO 100 

GLS02240 

29 CALL CALC { XP , YP , NSMALL , NB IG, P ) 

GLS02250 

DO 6 1=1, NSMALL 

GLS02260 

xi n=xpc i) 

GLS02270 

Y ( I )=YP( 1 ) 

GLS02280 

6 CONTINUE 

GLS02290 

RETURN 

GLS02300 

END 

GLS02310 
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noon 


$1 

c 


8FTC FLINFN LlST.OECK 



FL I 00000 

FUNCTION SUBROUTINE 

TO GENERATE INDEPENDENT VARIABLE 

FUNCTIONS 

FL I 000 10 

FUNCTION FUNC (I 

,X) 



FL I 00020 

z=x 




FL 100030 

K= I 




FL 100040 





FLI00050 

THE NUMBER UF THE STATEMENT IS IDENTICAL TO THE 

NUMBER 3F THE 

TERM FL I 00360 

IN THE POLYNUMIAL 

TO BE FITTED. 



FL 100070 





FL I 00080 

IF ( K-6 ) LOO, 6, 7 




FL 100090 

7 WRITE (6,10) K 




FLI00100 

10 F0RMAT(IH0,5X,72HT00 MANY COEFFICIENTS ASKED 

FOR 

- ONLY 6 

ARE AVFLI00110 

1AILABLE YOU CANNOT GET 13) 



FL 100120 

STOP 




FL I 001 30 

6 ANS=ALOG10 ( L ) 




FL 100140 

GO TO 99 




FL I 001 50 

100 IF (K-4) 101,4,5 




FL I 001 60 

5 ANS= Z*Z 




FLI00170 

ANS=ANS* ANS 




FL100180 

GO TO 99 




FLI00190 

4 ANS=Z*Z*Z 




FL I 00200 

GO TO 99 




FL I 002 10 

101 IF (K-2) 1,2,3 




FLI00220 

3 ANS =Z*Z 




FL 1 002 30 

GO TO 99 




FL 100240 

2 ANS=Z 




FL I 002 50 

GO TO 99 




FL 100260 

l ANS=1. 




FL 1002 70 

99 FUNC=ANS 




FL 100280 

RETURN 




FL 100290 

END 




FL 100300 
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r> n n n r> 


Y C I > ARRAY 


BL I 00300 
BL I 0001 0 
BL 1 00020 
BL 100030 
BLI00040 
BL 100050 

BE USED IN SPECIFYING BLI00060 

BLI00070 
BLI00080 
BL 100090 
BLI00100 
BL 1001 10 

10 FORMAT! 1H0 , 5X, 4 AH WRONG BACK-TRANSFORMATION OF Y REQUESTED, K= IJIBLI00120 


$1 BFTC BLINFN LIST, DECK 
C SUBROUTINE TO BACK-TRANSFORM 
FUNCTION BFUNC ( I , X) 

Z=X 
K= I 


THE NUMBER OF THE STATEMENT IS THE NUMBER TO 
THE TYPE OF CONVERSION. 

IT MUST AGREE WITH THE CONVERSION USED IN 

IF ( K-6 ) 100,6,7 
7 WRITE (6,101 K 


SUBROUTINE FUNC. 



STOP 

BL 1001 30 

6 

ANS=10.**Z 

BL 100140 


GO TO 99 

BL 1 001 50 

100 

IF ( K-4 ) 101,4,5 

BL I 00 l 60 

5 

ANS=SQRT ( Z ) 

BL I 001 70 


ANS=SQRT (ANSI 

BLI00180 


GO TO 99 

BLI00190 

4 

ANS=Z**. 333333333 

BL 100200 


GO TO 99 

BL 1002 10 

101 

IF ( K-2 ) 1,2,3 

BLI00220 

3 

ANS=SURT(Z) 

BLI00230 


GO TO 99 

BLI00240 

2 

ANS=Z 

BL 1 00250 


GO TO 99 

BL100260 

l 

ANS=l. 

BL 100270 

99 

BFUNC=ANS 

BLI00280 


RETURN 

BL 100290 


END 

BL 100300 
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BFTC DINV LIST, DECK, DEBUG 

OINOOOOO 

SUBR USING CROUT REDUCTION AND HENRY INVERSE ROUTINE 

DINOOOIO 


SUBROUTINE U INVER ( B, NNN , A , MAX, DE T , ERR, D ) 

D1N00320 


DIMENSION A I MAX, MAX) , B ( M AX , MAX ) , DC MAX ) , AN ( 60 ) 

DIN00030 


DOUBLE PRECISION A,B,D, AN 

DIN00060 


KM1=NNN 

DIN00050 


MM1=KM I— 1 

DIN00060 


HP1 = KMU1 

D1N00070 


DO 25 1=1, KM1 

DIN00080 


A ( I , MP 1 ) =DI I ) 

DIN00090 


DO 25 J= 1 , KM 1 

DIN00100 

25 

At I,J)=B(I,J) 

DIN00110 


DO 15 N=2,MP1 

DIN00120 

15 

At 1,N) =A( 1,N)/A( 1, 1) 

DIN00130 


DO 7 J=2 ,MP1 

DIN00140 


M = 0 

DIN00150 


L = J-l 

DIN00160 


DO 7 1=2, KM1 

DIN00L70 


M=M + 1 

D1N00180 


IF < M-L ) 10,10,9 

D1N00190 

9 

M = L 

DIN00200 

10 

DO 6 N= 1 , M 

DIN00210 

6 

At I,J)=At I,J)-A(N,J)*At 1,N) 

DIN00220 


IF(I-J) 12,7,7 

DIN00230 

12 

At I,J)=A( I , J ) / A ( 1,1) 

DIN00260 

7 

CONTINUE 

DIN00250 


DO 11 11=1, KM1 

DIN00260 


DEBUG (A(II,JJ),JJ=1,MP1) 

DIN00270 

11 

CONTINUE 

DIN00280 


0ET=1.0 

DIN00290 


DO 2 1 = 1, KM 1 

DIN00300 


Bt I,MPl)=At I, I ) *A t I,MP1)**2 

DIN00310 


DE T = DE T *A t 1,1) 

DIN00320 

2 

At I,I) = l./AtI,I) 

DIN00330 


DO 3J= 1 , MM 1 

DIN00340 


N=J*1 

DIN00350 


DO A I = N , KM 1 

DIN00360 


KJ = J 

DIN00370 


MJ= I-N 

DIN00380 


At I , J ) = A ( I,J)*A(J,J) 

DIN00390 


IF(MJ) A, A, 5 

D I NOO'tOO 

5 

DO 16 K= 1 , M J 

DIN00610 


KJ=K J*- 1 

D I NO 062 0 

16 

At I, J) =A( I , J)+At I ,KJ)*A< KJ, J) 

DIN00630 

A 

At I , J ) =- ( A ( 1 , J ) • A ( I, I) ) 

DIN00660 

3 

CONTINUE 

DIN00650 


DO 17 J = 1 , MM 1 

DIN00660 


KJ=KMl— J 

DIN00670 


DO 8 1 = 1, KJ 

DIN00680 


DO 8 K= 1 , J 

DIN00690 


KK=MP1-K 

DIN00500 

8 

At KJ, I ) =A(KJ, I ) - A ( KJ ,KK) *AtKK, I ) 

DIN00510 

17 

CONT INUE 

DIN00520 


DO 19 J= 1 , MM 1 

DIN00530 


KK=KMl— J 

DIN00560 


KJ=MPl- J 

DIN00550 


DO 110 I =K J , KM 1 

DIN00560 

no 

At MP1, I ) =A ( KK , 1 ) 

DIN00570 


DO 111 K =K J , KM l 

DIN00580 


At KK,K) =0. 

DIN00590 


DO 111 L=KJ , KM1 

D1N00600 
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111 A(KK,K)=A!KK,K)-A(MP1,L)*A!L,K) 

19 CONTINUE 

DO 88 11=1, KMl 

DEBUG! A ( II,JJ),JJ=1,MP1> 

88 CONTINUE 
REWIND 7 

WRITE (7) ( <B( I, J) , 1=1, KMl) , J=1,KM1) 
DO 114 I=1»KM1 
DO 112 J=1,KM1 
AN! J)=0. 

DO 112 K=1 , KMl 

112 AN! J)=AN!J)+B! I,K)»A(K,J) 

DO 114 1 1= l , KM 1 

B( I, 1 1 ) = (-AN( II ) ) 

IF ! 1 1- I I 114,113,114 

113 B( I, 1 1 ) =B! 1,11)4-1. 

114 CONTINUE 

DO 117 1=1, KMl 
DO 115 J=1,KM1 
AN! J)=0. 

UO 115 K=1 , KMl 

115 AN ( J ) = AN ( J ) +A ( I*K)*B!K*J) 

DO 116 11=1, KMl 

116 A! I , 1 1 ) = A! 1 , 1 1 ) +AN! II) 

117 CONTINUE 

DO 18 11=1, KMl 

DEBUG! A! II, JJ) , JJ = 1,MP1) 

DEBUG! B l II, JJ), JJ=1,MP1) 

18 CONTINUE 
REWINO 7 

READ 17) !!Bl I,J),I=1,KM1),J=1,KM1) 
ERR=0. 

DO 20 1=1, KMl 
DO 20 K= 1 , KMl 
ERR 1 = 0 • 

DO 21 J= 1 , KM l 

21 ERR1 = BI I,J)*A! J,K)*ERR1 
IF ! I— K ) 23,22,23 

22 ERR1=1 .-ERR1 

23 ERR1=ABS! ERR1 ) 

IF l ERR-ERR1 ) 24,20,20 

24 ERR=ERRi 

20 CONTINUE 
RETURN 
END 


DIN00610 
DIN00620 
DIN00630 
DIN00640 
DIN00650 
DIN00660 
DIN00670 
OIN00680 
DIN00690 
DIN00700 
DIN00710 
DIN00720 
DIN00730 
DIN00740 
OIN00750 
DIN00760 
DIN00770 
DIN00780 
DI NOO 790 
DIN00800 
DIN00810 
DIN00820 
DIN00830 
DIN00840 
DIN00850 
D I NOO 860 
DIN00870 
DIN00880 
DIN00890 
DIN00900 
DIN00910 
DIN00920 
0IN00930 
DIN00940 
DIN00950 
DIN00960 
DIN00970 
DIN00980 
DIN00990 
D I NO 1000 
DIN01010 
D I NO 1 020 
DIN01030 
DIN01040 
DIN01050 


21 



$1 BFTC DSL V L I ST , DECK , DEBUG 

C SUBROUTINE TO OBTAIN SOLUTION VECTOR FROM INVERSE 8 
C 

C SOLUTION VECTUR IS RETURNED IN THE AN ARRAY AND IN ((8(I,NP1) 
C WHICH IS THE AUuMENTEO MATRIX INVERTED 
C 

SUBROUTINE OMVMUL ( B, NNN , 0, AN , MA X ) 

DIMENSION B( MAX, MAXI ,D(MAX) ,ANI MAX) 

DOUBLE PRECISION 8,0, AN 

M=NNN 

NP1 = NNN-*1 

DO 100 1=1, M 

AN( I 1=0. 

DO 100 J = 1 , M 

100 AN( I)=B( I,J)*U(J)+AN( I) 

0011 11=1, M 

DEBUG! 8( II, JJ) , JJ=1,M) 

11 CONTINUE 

DEBUG! AN! I 1 1 , I 1 = 1, MI 
DO 4 1 = 1, NNN 
4 B( I , NP 1 ) =AN 1 I) 

WRITE 16,15) 

DO l A J= 1 , M 
DO 13 1=1, M 

13 B! MP 1 , I ) =B (I , J ) / SyRT ( B ( I,I)*B!J,J)) 

14 WRITE (6,16) ! B ( MP1 , K ) ,K=1,M) 

15 FORMAT (20H0C0RRELATI0N MATR I X= ) 

16 FORMAT! 1H08G16. 8) 

RETURN 

END 


DSLOOOOO 
OSLOODIO 
DSL00020 
1=1, NNN) OSLO 00 30 
DSL00040 
DSL00D50 
DSL00060 
DSL00070 
DSL00080 
DSL00090 
DSL00100 
DSLOOl 10 
DSL00120 
DSLOOl 30 
DSL00140 
DSLOOl 50 
DSL00160 
DSLOOl 70 
DSLOOl 80 
DSL00190 
DSL00200 
DSL00210 
DSL00220 
DSL00230 
DSL002 40 
DSL00250 
DSL00260 
DSL00270 
DSL00280 
DSL00290 
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non 


ilBFTC CAL LIST, DECK 

SUBROUT INfc CALC! X,Y,N,K,P) 

DETERMINE AN EMULSIUN CALIBRATION CURVE FOR INPUT TO LAGRANGE 

INTERPOLATION ROUTINE. 

COMMON IW 

DIMENSION XI250) , Y( 500 ) , P( 35) 

I XX= 2. 

XSTOP=95. 

X( 1 ) =XX 
DO 10 1=1,10 
N= I 

YY=P( 1 ) -fp( 2 ) *XX+P ( 3) *XX**2+P(^) *XX**3 +P(5)*XX*»4 
XX= YY 
XI 1+ 1 ) = Y Y 

IF (XX-XSTOP) 10,10,20 
10 CONTINUE 
20 K=N- 1 
NN=N 

4 DO 30 I = 1 , N 
Y ( I ) =2 . **NN 
NN=NN- 1 
30 CONTINUE 
N=N+1 
XI N ) =100 • 

YIN) =0. 


CAL00000 
CAL00010 
CAL00020 
CAL00030 
CAL00040 
CAL00050 
CAL00060 
CAL00070 
CAL00080 
C AL00090 
CAL00100 
CAL00110 
CAL00120 
CAL001 30 
CAL00140 
CAL00150 
CAL00160 
CAL00170 
CAL00180 
CAL00190 
CAL00200 
CAL00210 
CAL00220 
CAL00230 
CAL002<tO 
CAL00250 


WRITE! 6,99) CAL00260 

WR ITE ( 6, 103 ) (XII), Y( I),I = 1,N) CAL00270 

99 FORMAT! IHO, IOOH X VS. Y (DEFINES EMULSION CALIBRATION CAL00280 
XCURVE ) FOR INPUT TO LAGRANGE ROUTINE. ) CAL00290 


103 FORMAT! 1H0,2F10. A) 
IW=0 


CAL00300 

CAL00310 


101 RETURN 


CAL00320 


END 


CAL00330 
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$ I BFTC LGRN L I ST , DECK, DEBUG LGROOOOO 

SUBROUTINE LGRNGE ( X, Y, N, ARG, ANS, DER , KODE ) LGROOOIO 

DIMENSION X(200),Y(200),A(10),C(4> LGR00020 

EQUIVALENCE I XX, All) > , ( XO, A( 2 )) , { XI , A ( 3 ) > , ( X2 , A ( 4)) , ( X3 , A I 5 ) ) , LGR3UO40 
i(Y0,A<6l),(Yi,Ami,(Y2»A(8>l,(Y3 t A{9n,(CV,Cm>,{C2,Cm> , ( C3 , C ( LGR0005D 


10 

11 

13 

12 


15 

14 

51 

1 

2 

50 


23) ), l C4, C ( 4 ) ) 

1FIARG— X(2) ) 10 , 10* 11 
MX = 1 

GO TO 51 

IF ( ARG-X(N-l) ) 12,12,13 
MX=N- 3 
GO TO 51 
K = N- 1 

DO 14 JA=2,K 
IFIARG-X(JA) ) 15, 15,14 
MX=J A-2 
GO TO 51 
CONTINUE 
IF(KODE-l) 1,2, 1 
C0N=0 . 0 
GO TO 50 
CON=i. 

XX=ALOG( ARG+CON) 

DO 3 1=1,4 
MX I=MX+ I - 1 

A ( H-1)=AL0G(X(MXI)«-C0NI 
A( I*5)=AL00( Y< MX I l+CON) 

DEBUG XX,X0,X1,X2,X3 
ANS=i ( XX-X1 )*( XX-X2)*( XX-X3)/( ( XO- XI J * ( X0-X2 ) * ( X0-X3 ) ) )*Y3 


LGR00060 
LGR00070 
LGR30080 
LGR00090 
LGR00100 
LGR30110 
LGR00120 
LGR00130 
LGR00140 
LGR00150 
LGR001 60 
LGROO 170 
LGROO 180 
LGROO 1 90 
LGR00200 
LGR00210 
LGR00220 
LGR00230 
LGR00240 
LGR00250 
LGR00260 
LGR002 70 
LGR00280 
♦ t XX-X3LGR0Q290 


1)*(XX— X2)*(XX-X3) *Yl/( ( XI— X0)*( X1-X2) *( X1-X3) ) * ( XX-XO ) * ( XX-X1 I • LGR00300 
2 ( XX-X3 ) * Y2/ ( (X2-X0) »( X2-X1 )*{ X2-X3) )♦■( XX-XO)*( XX- XI) »1XX-X2)*Y3/ LGROO 3 10 
3 ( ( X3-X0 ) *( X3-X1) •{ X3-X2 ) ) LGROO 320 

ANS=EXP( ANSJ-CON LGR00330 

DEBUG ANS LGR00340 

GO TO (5,4), KQDE LGR00350 

4 DO 7 1=2,5 LGR00360 

7 C( I- 1 ) =XX-A ( I ) LGR00370 

CC0=Y0/(X0-Xl)/(X0-X2)/( X0-X3) LGROO 3 80 

CC1=Y1/ ( XI- XO) /( X1-X2) /( X1-X3) LGR00390 

CC2 = Y2/(X2-X0)/( X2-XU/I X2-X3) LGR00400 

CC3-Y3/(X3-X0)/(X3-Xl)/( X3-X2) LGR00410 

0ER=ANS/ARG*(CC0*(C2*C3+C2*C4+C3*C4)+CC1*( Cl*C3 + Ci*C4*C3»C4) «-CC2*( LGR00420 

1C1*C2+C1*C4+C2*C4)+CC3*( C1*C2+C1*C3+C2*C3) ) LGR00430 

5 RETURN LGR00440 

END LGR00450 
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noon r r> o o o 


$ I B F TC RAT LIST, DECK RATOOOOO 

SUBROUTINE fATIO (NEL,Z,ZB) RAT00010 

RAT 00020 

COMPUTES RATIO UF KcLATIVE INTENSITY OF AN ELEMENT TO RELA1IVE RAT00030 
INTENSITY OF A STANDARD, MAXIMUM NUMBER OF STANDARDS IS FOUR, RAT00040 

MAXIMUM NUMBfR OF ELEMENTS IS FIFTY. RAT00050 

RAT 0C060 

REAL MIS RAT00070 

COMMON I w , KS T OR , NEXF RAT00080 

COMMON Z( 50, 1 I, 7) , ZB( 50, LI, 7) RAT00090 

COMMON EL ( 50 ) RATOOIOO 

DIMENSION RAT (. 1 1 , 4 , 7 ) , XR A T{1 1 , 4 , 7 1 , A VG ( i I , 4 ) , X A VG ( 1 1 , A ) RAT00110 

DATA R IEN, ETL ILE/6H , 3H* / RATC1200 

DIMENSION MIS( 11, A, 7) RAT00130 

WR I TE ( 6 , 98 ) RA T 00 140 

WRITE! 6,99) ( TLI I ), 1 = 1 ,KSTOP I RAT00150 

WRITE! 6, 101) RAT 00 160 

WRITE (6, 105) RAT 001 70 

WRITE16.106) RAT00180 

KST=K5T0P*1 RAT00190 

DO 40 1=KST , NEL RAT00200 

DO 30 KK= 1 , K STOP RAT00210 

DO 20 J = 1 , N E X ? RAT 00 220 

M=I RAT00230 

XN = C • RAT00240 

XNN =0 . RAT00250 

XAVG! J » KK ) =0 • RAT 00260 

A VG! J » KK ) =C . RAT00270 

SUM=0. RAT00280 

SUM 1 = 0 • RAT00290 

DO 10 K = 1,7 RAT00300 

N = K RAT 00 31 0 

L=K RAT 00 320 

RAT! J,KK,K)=0. RAT00330 

XR AT ( J , KK, K ) =0 . RAT00340 

IF! I-KST ) 1 » 4 , 1 
4 MIS! J,KK,K)=RIEN 


ARRAY MIS SHOWS PROGRAM EXTRAPOLATED A REL INTENSITY 
NO ORIGINAL DATA FOR THAT EXPOSURE AND STEP. 

1 IF (Z(I,J,K) >2,6,2 

6 IF (MIS! J,KK,K1-ET0ILE) 10,7,10 

7 Z(KK,J,K)=0. 

MIS! J,KK,K)=R IEN 
GO TO 10 


IF THERE 


RAT 00360 
WASRAT00370 
RAT 00380 
RAT00390 



2 IF ( Z ( KK, J , K » ) 3,11,3 

RAT 00410 

c 

TEST FOR BACKGROUND 

RAT 00420 


3 IF ( ZB ( I , J , K ) I 32,5,32 

RAT 00430 


5 IF ( ZB ( KK , J , K ) I 32,31,32 

RAT 00440 

c 

IF IN STEP 7 GU TO 21 AND START SEARCHING LOWER STEPS. 

RAT0D450 

c 

IF NOT IN STEP 7 START SEARCHING HIGHER STEPS. 

RAT00460 


11 IF (L-7) 12,21,21 

RAT00470 


12 L = L+l 

RAT 00480 


N = L-K 

RAT 00490 


IF(Z(KK,J,L) ) 13,11,13 

RAT00500 

c 

CORRECT FOR STEP MOVEMENTS 

RAT00510 


13 Z ( K K , J , K ) = Z ( KK , J , L ) * 2.**N 

RAT00520 


MIS! J,KK,K)=F IOILE 

RAT00530 


L4 IF ( ZB ( KK , J , L ) ) 15,3, 15 

RAT00540 


15 Z B ( KK , J , K ) = 7iJ(KK,J,U • 2.**N 

RAT00550 


GO TO 32 

RAT 00560 

c 

START SEARCH UF LOWER STEPS 

RAT 00570 
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non 


21 

22 


31 

33 
10 

35 

34 

36 
20 
30 


47 


46 


45 


- ZB! I, J,K) )/(ZIKK, J,K) - Zt)< KK, J,K) ) 
/ Z(KK,J,K> 


45 


J=1,NEXP,KLJ 


I F ( L— 7 ) 11,22,22 

L = L-l 

I F ( Z ( KK , J , L ) ) 23,27,23 
23 N = K-L 

CORRECT FOR STEP OUVEMEnT 

Z(KK,J,K)=Z(KK,J,L)/2.*»N 
MIS! J,KK,K)=EtOILE 
ZB(KK,J*K)=2b<KK,J,L>/2.**N 
I F t Z8(KK,J,L)> 32,3,32 
27 IF (L-l) 10,10,22 
DETERMINE RATIO 
32 RAT ( J,KK,K) = ( Z< I, J,K) 

XN=Xn+1. 

XRAT ( J , KK, K ) = Z ( I , , K ) 

XN.N = XNN + 1. 

SUM = SUM «■ RaT(J,KK,K) 

SUM1=SUM1+XRA1 ( J,KK,K) 

CUNT INUE 
I F ( SUM) 35, 34, 35 
AVGtJ.KK) = SUM/ XN 
IF! SUM 1 ) 36,20,36 
XAVG 1 J , KK ) = SUM1/XNN 
CONI INUE 
CUN T I Nut 

WRITE (6, 102) ( ELI I ) ,EL(KK) ,KK=1 ,KSTOP) 

IFIKSTUP-2 )47 , 47 ,46 
KLJ=2 
NN = 0 
GO TU 
KLJ = 1 
NN= 1 
DO 50 
LJ=J+1 

WRITLt 6, 104) 

DO 51 J J = 1 , 7 
WRITEI6,100) 

XKSTOP ) 

1FINN) 51,53,51 

53 WRITE(6,1C7) L J , J J , ( RAT ( L J , KK , J J ) , XRAT ( L J , KK , J J ) ,MIStLJ,KK,JJ 
X= 1 » KS TOP ) 

51 CONTINUE 

WR HE ( 6, 103) ( AVGtJ.KK) ,XAVC(J,KK) ,KK = l,KSTOP) 

IFINN) 50,48,50 

48 WR I T E ( 6 , 108 ) ( A VG t L J , KK ) , XA VG ( L J , KK ) , KK = 1 , K STOP ) 

50 CONTINUE 
40 CONTINUE 

FORMAT STATEMENTS 

98 FORMAT (lHl) 

99 FORMAT! 1H2.57X, 16HINTENSITY RATIOS //.37X.59H THE FOLLOWING. 

XNTS Are CONSIDERED INTtRNAL STANDARDS- , / / , 46 X , ( 4 ( A6 , 4X )) ) 

100 FORMAT t 1H , 2 1 4 , 4X , t 4 ( 2G 1 1 . 4, A 3 ))) 

101 FORMAT! 1H0, 15X, 105H COLUMNS 1,4,7,10 ARE BACKGROUND CORRECTED 

XOS. COLS. 2,5,8,11 ARE THE UNCORRECTFD RATIOS. ) 

102 FORMAT! 1H1, 12H EXP STEP ,X,4IA6,2h /,A6,11X)) 

103 FORMAT ( 1H , 5X, 3MAVG.4X, ( 4( 2G1 1 .4, 3XJ) ) 

104 FORMAT (1H0) 

105 FORMAT l 1H0, 35X , 60H * IN COL. 3,6,9,12 MEANS VALUE OF STD. IS 

XPOLATED. ) 

106 FORMAT! IHO , 32X, 64H A BLANK IN THESE COLUMNS MEANS RATIO WAS 

XLATED DIRECTLY. ) 

107 FORMAT! 1H+60X, 2 14, < 2{ 2G1 1 . 4, A3 ) ) ) 

108 FORMAT ( 1H+64X, 3HAVG, 2X, ( 21 2GL I .4, 3 X ) ) ) 

RETURN 

END 


J,JJ, (RAT 1 J , KK , J J ) , XR A T ! J , KK , J J ) , M I S ( J , KK , J J ) , KK= 1 , 


RAT00580 
RAT 00590 
RAT00600 
RAT00610 
RAT 00620 
RAT 00630 
RAT00640 
RAT 00650 
RAT00660 
RAT 00670 
RAT00680 
RAT 00690 
RAT00700 
RAT00710 
RAT00720 
RAT00730 
RAT00740 
RAT00750 
RAT 00760 
RAT 00770 
RAT 00780 
RAT00790 
RAT00800 
RAT 00810 
RAT00820 
RAT00830 
R AT00840 
RAT00850 
RAT 00860 
RAT00870 
RAT00880 
RAT 00890 
RAT00900 
RAT 0091 0 
RAT 00920 
RAT 00930 
RAT0094C 
RAT00950 
) , KK RAT0096C' 
RAT 00970 
RAT 00980 
RAT00990 
RAT01000 
RAT 01010 
RATO 1020 
RAT01030 
RAT 01040 
RATO 1050 
RAT 01060 
RAT01070 
ELEMERAI09300 
RAT09350 
RAT 09380 
RATI RAT 01110 
RAT01120 
RAT09500 
RATO 1 1 40 
RATO 1 1 50 
EXTRARAT01 160 
RATO 1170 
CALCURAT01 180 
RAT01 190 
RATO 1 200 

RATO 1 220 
RA 1331230 
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APPENDIX C 


SAMPLE OUTPUT FROM SUBROUTINE RATIO 
The following output is a typical example from the subroutine RATIO: 

INTENSITY RATIOS 

THE FOLLOWING ELEMENTS ARE CONSIDERED INTERNAL STANDARDS- 
A 40AA W A008 AG2929 AG3501 

COLUMNS 1,9,7,10 ARE BACKGROUND CORRECTED RATIOS. COLS. 2,5,8,11 ARE THE UNCORRECTED RATIOS. 
* IN COL. 3,6,9,12 MEANS VALUE UF STD. IS EXTRAPOLATED. 

A BLANK IN THESE COLUMNS MEANS RATIO HAS CALCULATED DIRECTLY. 


EXP 

STEP 

N I 3566 

/A 4044 


NI3566 /W 4008 

N 13566 f AG2929 

N 1 3566 

/AG3501 


1 

1 

C.1459 

C. 1743 

* 

32.655 5.0638 

0.3581 0.4277 

0.1020 

0.1218 

* 

1 

2 

C 

0 


0 0 

0 0 

0 

0 


i 

3 

G 

c 


0 0 

0 0 

0 

0 


1 

4 

0 

0 


0 0 

0 0 

0 

0 


1 

5 

C 

0 


0 0 

0 0 

0 

0 


1 

6 

c 

0 


0 0 

0 0 

0 

0 


1 

7 

c 

0 


0 0 

0 0 

0 

0 



AVG 

0.145S 

0. 1743 


32.655 5.0638 

0.3581 0.4277 

0.1020 

0.1218 



2 

1 

C 

c 


0 

0 

0 

0 


0 

0 


2 

2 

0 

0.8380 

♦ 

0 

0.8663E-01* 

0 

0.2285 

* 

0 

0.1537 

* 

2 

3 

C 

C. 8498 


0 

0 • 8785 E— 01* 

0 

0.2317 

* 

0 

0. 1558 

* 

2 

4 

C 

0.9295 


0 

0.8646E-01* 

0 

0.2281 

* 

0 

0.1534 

♦ 

2 

5 

c 

0 


0 

0 

0 

0 


0 

0 


2 

6 

0 

0 


0 

0 

0 

0 


0 

0 


2 

7 

c 

0 


0 

0 

0 

0 


0 

0 



AVG 

c 

0.8724 


0 

0.8698E-01 

0 

0.2294 


0 

0.1543 



3 

l 

C 

c 


0 

0 


0 

0 

0 

0 

3 

2 

C 

0 


0 

0 


0 

0 

0 

0 

3 

3 

C 

0 


0 

0 


0 

0 

0 

0 

3 

4 

c 

2-9449 

* 

0 

0.1648 

* 

0 

0.8834 

0 

0. 6593 

3 

5 

0 

3.1737 

* 

0 

0.1776 

* 

0 

0.9145 

0 

0. 7719 

3 

6 

0 

0 


0 

0 


0 

0 

0 

0 

3 

7 

c 

0 


0 

0 


0 

0 

0 

0 


AVG 

0 

3. 0593 


0 

0.1712 


0 

0.8990 

0 

0.7156 


4 

1 

C.2297 

0.2774 

♦ 

0.3154E-01 

0.3810E— 01* 

0.6918E-01 

0.8356E-01* 

0.589 4E-01 

0.711 9E-01 * 

4 

2 

0.2630 

0.2818 


0. 361 IE— 01 

0 . 3870E-0 1* 

0.792QE-01 

0.8488E-01* 

0.6748E-01 

0. 7231E-01* 

4 

3 

C 

0 


0 

0 

0 

0 

0 

0 

4 

4 

G 

0 


0 

0 

0 

0 

0 

0 

4 

5 

C 

0 


0 

0 

0 

0 

0 

0 

4 

6 

0 

0 


0 

0 

0 

0 

0 

0 

4 

7 

c 

c 


0 

0 

0 

0 

0 

0 


AVG 

0.2463 

0.2796 


0.3383E-01 

0.3840E-01 

Q.7419E-01 

0.8422E-0 1 

0.632 IE-01 

0.71 75E-01 
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5 

1 

C.5671E-C1 

0.1359 

* 

0.7482E-G2 

0.1052E-01* 

0.4045E-01 

0.5685E-01* 

0.3139E-01 

0*44126—01 * 

5 

2 

0. 1244 

0. 1379 


0.9628E-02 

0.1067E-01* 

0.5206E-01 

0.5769E-01* 

0.4040E-01 

0.4477E-01* 

5 

3 

C 

0 


0 

0 

0 

0 

0 

0 

5 

4 

0 

G 


0 

0 

0 

0 

0 

0 

5 

5 

G 

C 


0 

0 

0 

0 

0 

0 

5 

6 

0 

0 


0 

0 

0 

0 

0 

0 

5 

7 

c 

C 


0 

0 

0 

0 

0 

0 


AVG 

C • 1 1 06 

0. 1369 


0.8555E-G2 

0 . 1 059E-01 

0.4625E-01 

0.5727E— Oi 

0. 3590E-01 

0.4445E-01 


6 

1 

C 

c 


0 

0 


0 

0 

0 

0 

6 

2 

C 

0 


0 

0 


0 

0 

0 

0 

6 

3 

0 

0 


0 

0 


0 

0 

0 

0 

6 

4 

c 

4.9703 

♦ 

0 

0.7166 

• 

0 

0.9195 

0 

0.6623 

6 

b 

0 

5.0939 

* 

0 

0. 7344 


0 

0.8736 

0 

0. 7682 

6 

6 

0 

0 


0 

0 


0 

0 

0 

0 

6 

7 

0 

0 


0 

0 


0 

0 

0 

0 


AVG 

c 

5.0321 


0 

0.7255 


0 

0.8965 

0 

0.7153 


7 

1 

0 

0 

0 

0 


0 

0 


0 

0 


7 

2 

0 

0.8036 

0 

0.1225 

* 

0 

0.2686 

* 

0 

0.2009 

* 

7 

3 

c 

1.0517 

0 

0.1363 

♦ 

0 

0.2988 

* 

0 

0.2235 

* 

7 

4 

0 

0 

0 

0 


0 

0 


0 

0 


7 

5 

c 

0 

0 

0 


0 

0 


0 

0 


7 

6 

c 

0 

0 

0 


0 

0 


0 

0 


7 

7 

0 

0 

0 

0 


0 

0 


0 

0 



AVG 

0 

C.9276 

0 

0. 1294 


0 

0.2837 


0 

0.2122 



8 

1 

C 

C.4C8C 

* 

0 

0.9922E-01* 

0 

0.1167 

* 

0 

C. 9206E-0 1 * 

8 

2 

C 

0.4189 

* 

0 

0.10L9 * 

0 

0.1199 

* 

0 

0.9452E-01* 

8 

3 

c 

0 


0 

0 

0 

0 


0 

0 

8 

4 

c 

0 


0 

0 

0 

0 


0 

0 

8 

5 

c 

c 


0 

0 

0 

0 


0 

0 

8 

6 

c 

0 


0 

0 

0 

0 


0 

0 

8 

7 

0 

c 


0 

0 

0 

0 


0 

0 


AVG 

0 

0.4135 


0 

0.1005 

0 

0.1183 


0 

0.9329E-01 
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